Mapping of Coulomb gases and sine-Gordon models to statistics of random surfaces 
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We introduce a new class of sine-Gordon models, for which the interaction term is present in a 
region different from the domain over which the quadratic part is defined. We develop a nonpertur- 
bative approach for calculating partition functions of such models, which relies on mapping them to 
statistical properties of random surfaces. As a specific application of our method, we consider the 
problem of calculating the amplitude of interference fringes in experiments with two independent 
low dimensional Bose gases. We calculate full distribution functions of interference amplitude for 
one-dimensional and two-dimensional gases with nonzero temperatures. 



I. INTRODUCTION 

Sine-Gordon (SG) models and their Coulomb gas rep- 
resentations appear as a low energy effective theory for 
many types of physical systems. The so-called "bulk" 
SG model describes the Berezinskii-Kosterlitz-Thouless 
(BKT) transition in two-dimensional superfluids [l|, 0, HJ 
and the superfiuid to insulator transition of Cooper 
pairs in a chain of Josephson junctions 4]. The so- 
called "boundary" SG model @ can be used to describe 
the Chakravarty-Schmid transition in a single Josephson 
junction with dissipation [6| and a quantum impurity 
problem Q in one dimension (ID). Powerful theoreti- 
cal techniques have been developed for studying such SG 
models including Bethe ansatz solutions [8(, renormal- 
ization group analysis (see e.g., Ref. [3[), and functional 
renormalization group [9(. Recent theoretical work sug- 
gested that there is another class of SG models which 
is important for several kinds of physical systems. Such 
models can be described by the action 

S(g)=irK I (W4>) 2 dxdT + 2 j g cos (27rcp)dxdr, (1) 

where the interaction term cos (2n(j)) is present in the 
spatial region ui (or space-time region for quantum prob- 
lems), which is only a part of the domain ft over which 
the noninteracting part (V(/>) 2 is defined. This should be 
contrasted to the bulk SG models, in which the interac- 
tion term is present in the entire region ft, and to the 
boundary SG models, in which the interaction term is 
present on a line. We note that models with inhomoge- 
neous g(x, t) can be considered using methods developed 
in the paper as well, but for concreteness we will consider 
only constant g. Model |T]) interpolates between the bulk 
and boundary SG models, and we will refer to it as the 
interior sine-Gordon (ISG) model. Here are a few exam- 
ples of physical systems that can be described by such 
models. The first example is the problem of "interwire 
coherence" [To| in which wires are brought together over a 
finite length I and separated on both ends. The cos-term 
describes the correlated umklapp electron scattering in 
the two wires (+2kp scattering in one wire accompanied 
by —2kp scattering in the other wire). The quantum 



space-time action of this system has the form of Eq. (fTJ) 
with ft = [-oo,oo] x x[0,$ T andw = [-1/2, l/2] x x [0, [3} T . 
Another example comes from a system of quantum par- 
ticles in one dimension (e.g., electrons or Cooper pairs in 
a wire, or ultracold atoms in a weak optical trap) with a 
periodic potential present in a finite region of the system. 
The cos-term comes from the umklapp scattering on the 
external potential and is limited to the finite region in the 
interior of the system. The third example is the problem 
of distribution functions of interference fringe amplitudes 
(DFIFA) for apair of independent low-dimensional con- 
densates [ll], Ell [H| ■ Individual moments of the distri- 
bution function can be represented as a microcanonical 
partition function of Coulomb gases 0, [TH, fig , with 
positions of Coulomb charges restricted to the part of 
the system from which interference patterns have been 
extracted. The latter is typically smaller than the total 
system size. For example, in the case of large 2D con- 
densates we have ft = [—00,00] a x [— 00,00] T , and when 
the interference pattern is extracted from the area Z x Z, 
we have w = [— l/2,l/2] x x [— Z/2,Z/2] T . The relation be- 
tween the DFIFA and the partition function of the ISG 
model will be outlined below and has also been discussed 
inRefs. dlEEm. 

In this paper we develop a nonperturbative approach 
to calculating partition functions of a wide class of SG 
models and Coulomb gases, which relies on the mapping 
of their partition functions to certain problems of statis- 
tics of random surfaces. We point out that our method 
does not rely on the existence of the exact solutions of 
SG models, but uses the structure of the multi-point cor- 
relation functions in the absence of interactions. We also 
note that a suitable extension of our method can be used 
to compute correlation functions of SG models in equi- 
librium and non-equilibrium situations. The particular 
strength of our approach is that it can be applied to study 
ISG models, which cannot be analyzed by other theoreti- 
cal methods. As a concrete application of our method we 
calculate DFIFA for both ID and 2D condensates. We 
point out that earlier theoretical work on interference ex- 
periments focused on ID systems with periodic boundary 
conditions (PBC) [l5| . While these boundary conditions 
are extremely artificial from the point of view of realis- 
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tic experiments, they allow one to relate the DFIFA to 
the quantum impurity problem [?J and use certain exact 
results about the latter [T3|. Methods used in Ref. [l5[ 
can not be generalized either to the realistic case of open 
boundary conditions (OBC) (i.e., interference patterns 
extracted from the interior of a large system) or 2D sys- 
tems, but these cases can be analyzed using the method 
discussed in this paper. We emphasize, however, that the 
main goal of this paper is to introduce an approach to 
the analysis of SG models and the problem of DFIFA is 
just one example that illustrates the power of the new 



method. 



II. MAPPING 

The partition function corresponding to the action (TTJ) 
is given by Z(g) = jV^e~ s ^/ /2?0e- s(o) . By expand- 
ing Z(g) in powers of g we arrive at the grand canonical 
partition function of the Coulomb gas [18| 



Z(g) = -hvl 2 ^ where Z *n = ■■■ d 2 u 1 ...d 2 v n e V<* « > . (2) 
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Here, Zi n is a microcanonical partition function of a 
classical two-component neutral Coloumb gas of In par- 
ticles, and G(x, y) is an interaction potential, which is 
proportional to Green's function of the Laplace operator 
on Q. The most familiar case is when fi = [— oo,oo]a; x 
[—00, oo] r and G(x, y) = ln\x — y\. 

To evaluate Z(g) nonperturbatively in g, we introduce 
an auxiliary function W{oi), which should be understood 
as a certain distribution function, and is defined in such 
a way that its nth moment equals Z 2 „. 



Z 2n = / W(a)a n da. 



(3) 



One can use the Hankel transformation [19( to compute 
Z(g) from VF(oi) as 



z(g) 



W(a)I (2g^)dc 



(4) 



This equation equation can be verified using the Taylor 
expansion of the modified Bessel function I (x) and Eq. 
© . Formulation of the auxiliary " problem of moments " 
allows one to avoid calculating Z(g) order by order, and 
can be viewed as a tool to sum the perturbation series in 
Eq. © to all orders. 

Function G(x, y) is real and symmetric, so it can be 
diagonalized on to by solving the eigenvalue equations 



G(x,y)^ f (y)d 2 y = G(f)^ f (x). 



(5) 



Here / is an integer index, which goes from 1 to 00. ^/(x) 
can be chosen to be real and normalized according to 
J^^ f (x)^ k {x)d 2 x = 5(f,k). Then, G(x,y) is given by 

G(x,y) = J2f=T Such decomposition 

is similar to the diagonalization of a symmetric matrix 
using its eigenvectors and eigenvalues. We have 



d 2 u l ...d 2 v n e Y ' ! ^ */( s -)-*/(^)) 2 -x:'=r(*/(*) 2 +*/(^ 



d 2 ui...d 2 v n Yl 



J —OO J 



2tt 



(6) 



To go from the first to the second line in Eq. (J6j) we u— and v— integrals are identical, 
introduced the Hubbard- Stratonovich variables tf. Inte- 
gration over d 2 u\...d 2 v n , is now straightforward since all 




9({tf}) n g({-t f }r, (7) 
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where 

g({t f }) = [ dx e^/'/V^F 1 */^-^*/^ 2 . (8) 

From a comparison of Eqs. ([3]) and {7J we obtain 

f=oo roo e -% dt 

W(a)= J] J -°° f S[a-g({t f })g({-t f })]. (9) 

/=l v^vr 

Equations ([8|) and ([9]) have a simple physical interpre- 
tation. Consider ^Sf(x) to be the eigenmodes of the 
surface vibrations, tf the fluctuating mode amplitudes, 
and |G(/)| the noise power. Infinite dimensional inte- 
gral over {tf} variables can be understood as an aver- 
aging over fluctuations of the surface. For a particu- 
lar realization of noise variables {tf}, a complex valued 
surface coordinate at point x is given by h(x;{tf}) — 
Y, f tf^f{x)^G{f)/K-G{f)^f{x) 2 /2K. For each real- 
ization of a random surface {tf}, g{{tf}) is obtained as 
an integral (|8|), which can also be written as g({tf}) = 
d 2 xe h( - s '^ tf ^ . Hence Eq. jH) can be interpreted as the 
mapping between the partition function of the SG model 
([1]) and the statistics of random surfaces subject to clas- 
sical noise. This mapping is the central result of this 
paper. 

In general, g({tf}) is a complex number, thus in gen- 
eral, a is defined on a complex plane. Simplifications oc- 
cur: if all eigenvalues G(f) are negative, then tf}) = 
9{{tf})*,g{{tf})g({-t f }) = \g{{t f })\ 2 , and a is always 
real and positive. Function W(a) can be computed ef- 
ficiently using Eq. ((9]) and Monte Carlo simulations. In 
the first step one solves integral equations ([5|) numer- 
ically to obtain eigcnfunctions ^f(x) and eigenvectors 
G(f). Then one samples random numbers {tf} from the 
Gaussian ensemble, and plots the histogram of the results 
for g({tf})g({— tf}). Each point on a histogram requires 
a computation of only two integrals, and IT (a) can be 
evaluated to arbitrary precision. 

Note, that this simple numerical evaluation of PF(a) 
allows one to extract Z(g) for all values of g using Eq. 
((3]). While conventional large-scale Monte Carlo simu- 
lations [H HH can be used to extract properties of SG 
models, they require separate simulation for each value 
of g. Application of such methods to the calculation of 
DFIFA described below would also require analytic con- 
tinuation of numerical results to imaginary values of g 
[l5j . which is numerically unstable. In addition, the map- 
ping in Eqs. ([5]) and @ does not only simplify numerical 
simulations, but can also be used to obtain analytical 
results (see below). 

III. APPLICATIONS TO INTERFERENCE OF 
LOW-DIMENSIONAL GASES 

We now apply a general formalism developed to a par- 
ticular problem of the interference of low-dimensional 
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FIG. 1: Simplified setup of interference experiments with ID 
Bose liquids (see, e.g., Refs. [TTI. H3p. Two parallel conden- 
sates are extended in the x direction. After atoms are released 
from the trap, clouds are imaged by the laser beam propagat- 
ing along the z axis. Meandering structure of the interference 
pattern arises from phase fluctuations along the condensates. 
The net interference amplitude A is defined from the density 
integrated along the section of length L. For the analogous 
2D setup, see, e.g., Fig. 2 of Ref. [Q. 

Bose gases [H El El El EE 13 • Typical experimental 
setup for interference of low-dimensional gases with open 
boundary conditions (OBC) is shown in Fig. [TJ Two par- 
allel condensates are extended in the x direction. After 
atoms are released from the trap, clouds expand predom- 
inantly in the transverse direction. After sufficient time 
of flight clouds overlap, and the laser beam propagating 
along the z axis takes an absorption image. Fluctuations 
of the relative phase result in fluctuations of the minima 
positions for different x. For each y, the image can be 
integrated along the x direction to obtain the integrated 
fringe amplitude A. One experimental image can be used 
to extract information for different values of L. Many im- 
ages are still required to obtain distribution functions for 
each L. 

For two identical ID clouds higher moments of the 
fringe amplitude A can be written as [Tij 

(\A\ 2n ) = A 2n Z 2n , where A = ^Cp 2 ^^ L 2 ^/ K . (10) 

Here p is the density, L is the imaging length, £/j is 
the healing length, and C is a constant of the order 
of unity. For OBC Z 2n is given [ll| by Eq. (gj) with 
fi = [— oo,oo]j; x [0,/3] r (and periodic boundary con- 
ditions in t) and uj — [0,1]^ (and fixed r). Equation 
(HHJ) has been derived neglecting the shot noise, which 
arises due to a finite number of particles in the inter- 
fering clouds 0, [13]. In what follows, we will be in- 
terested in the distribution functions W(a) of a pos- 
itive variable a — \A\ 2 /Aq defined by Eq. ([3]), or 
of its normalized version a = \A\ 2 /(A 2 ), defined by 
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FIG. 2: Distribution functions of the normalized interfer- 
ence amplitude W(a) at T = for ID gases with open 
boundary conditions, shown for Luttinger parameters K = 
2 (dashed), if = 3 (dotted), and K = 5 (solid). The in- 
set shows a comparison between open (solid) and periodic 
(dashed) boundary conditions for K — 5. 
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FIG. 3: Distribution functions of the normalized interference 
amplitude W(a) for a ID Bose gas with open boundary condi- 
tions at nonzero temperature and K = 5. Different curves cor- 
respond to ratios K^t/L — oo (solid) , K^t/L = 1 (dotted), 
and K^t/L = 0.25 (dashed). £t is the thermal correlation 
length, K is the Luttinger parameter, and L is the imaging 
length. 



Z 2n /^2 = /o°° W(a)a n da. For zero temperature Z 2n 
depends only on the Luttinger parameter K, which de- 
scribes [23|, H|| the long-distance behavior of boson cor- 
relation functions, given by (a^(x)a(O)) ~ p(S,h/x) 1 ^ 2K . 
For bosons K ranges from K = 1 (strong interactions) to 
K = 00 (weak interactions). For OBC and zero temper- 
ature G(x, y) equals G(x, y) = log \x — y \ , while for PBC 
considered in Ref. 15} , G per (x,y) = In — sin7rjx — y\. 
While A depends on L, for zero temperature distribu- 
tion W(a) does not depend on L. but depends only on 
K. For nonzero temperature, Zm depends on K and the 
thermal length £t = hv s /(kBT), where v s is the sound 

velocity: G{x, y, £ T /L) = In ( sinh 7r|a jg /|£ ) . 

For 2D, one can use a similar approach to describe 
the contrast distribution at finite temperature below 
the BKT transition. In this case, correlation functions 
are given by (a t (r)a(0)) ~ p(£h/r) , where rj(T) — 
mT / (2nh 2 p s (T)) depends on the temperature and the 
superfluid density p s (T). The BKT transition happens 
at the universal value rj c (T c ) = 1/4. To keep a connec- 
tion to the ID case, we will use K = l/(2rj(T)), and 
restrict our attention to K > K c = 2. For 2D with the 
aspect ratio of the imaging area equal to unity, Ui and 
Hi in Eq. @ are defined on a square ui = [0, l] x x [0, 1] T 
with G(x, y) = In \S — y\. 

In the limit K — > 00, one can expand the exponent of 
Eq. © in the Taylor series. Then a = g({tf})g({-tf}) 
is linearly related to roughness, or mea n sq uare fluctu- 
ation of the surface, as defined in Ref. [25|]. For PBC 
the noise has a 1/f power spectrum, since G per (x,y) — 
In i sin7r|x — y\ = — \n2n — J2f=T (cos2nfxcos27rfy + 
sin 2irfx sin 2ivfy) / f. This results in the Gumbel statis- 
tics of the roughness [H GJ, H S3]: W G (a) = 
K exp(x — e x ), where x — K{a — 1) — 7 and 7 = 0.577 
is the Euler constant. This provides the analytical proof 
of the conjecture made in Ref. [l5j], that the distribution 
function in this case is given by one of the extreme value 
statistical distributions, the Gumbel function [26| . 

In what follows we perform simulations of VF(a) with 



up to N — 10 6 — 10 7 realizations of {tf} and smoothen 
the data. We use a finite value of f max and check for 
convergence with f m ax, typically ~ 30. (a) is always kept 
within 1% from its expected value. For most of the pre- 
sented results, all eigenvalues G(f) are negative, and Eq. 
© can be directly applied. Special care should be taken 
of the ID case with nonzero temperatures, since one of 
the eigenvalues can be positive. This situation can be 
handled by subtracting a sufficiently large positive con- 
stant C from G{x,y,^T/ L), which makes all eigenvalues 
negative. According to Eqs. ([2]) and ([3]), this leads to 
rescaling of a by a factor e~°l K , which can be easily 
taken into account. 

In Fig. [2] we show distribution functions of the normal- 
ized interference amplitude W(a) at T = for ID gases 
with OBC for various K. The inset shows a comparison 
between OBC and PBC for K = 5. In Fig. El we show 
distribution functions of the normalized interference am- 
plitude for a ID gas with OBC at nonzero temperature 
and K = 5. For ^tK/L <C 1 distribution is Poissonian 
0GJ,Bi| and wide, while for K > 1 and £, T K/L > 1 it 
is very narrow. Evolution of the full distribution function 
of the visibilities as L is varied can be used to precisely 
measure the thermal length and to extract the tem- 
perature. As seen in Fig. (3J at T 7^ the distribution 
function has characteristic features, i.e., it is generally 
nonsymmetric and can have a minimum. These features 
can be used to distinguish the intrinsic noise due to fluc- 
tuations of the phase from technical noise. Finally, in 
Fig. [5] we show distribution functions of the normalized 
interference amplitude for a 2D gas with an aspect ra- 
tio of imaging area equal to unity and OBC below the 
BKT temperature. Above the BKT temperature, distri- 
bution functions become Poissonian for L > ^ where £ 
is the correlation length. In 2D one cannot describe the 
crossover at L ~ £ similar to ID, since the action which 
describes the fluctuations of the phase is not quadratic 
in this region, and Eq. © does not hold. 
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FIG. 4: Distribution functions of the normalized interfer- 
ence amplitude VF(<5) for a two-dimensional Bose gas with 
the aspect ratio of the imaging area equal to unity and open 
boundary conditions. Temperature is below the Berezinskii- 
Kosterlitz-Thouless (BKT) transition temperature. Different 
curves correspond to r}(T) — rj c (T c ) — l/4(the BKT transi- 
tion point, solid), 77 (T) = 1/6 (dashed line), and 77 (T) = 1/10 
(dotted line). Above the BKT transition temperature the 
distribution function is Poissonian (dot-dashed line). 



tial region different from the domain over which the non- 
interacting part is defined. We developed a general map- 
ping of such sine-Gordon models and related Coulomb 
gases to statistical properties of random surfaces, which 
can be used to calculate their partition functions non- 
perturbatively. As a specific application of our approach, 
we considered interference experiments with two inde- 
pendent low-dimensional Bose gases. We calculated full 
distribution functions of the interference amplitude for 
ID and 2D gases with open boundary conditions and 
nonzero temperatures. Full distribution functions of in- 
terference fringe visibilities can be used for thermometry. 



IV. CONCLUSIONS 

To summarize, we introduced a class of sine-Gordon 
models, for which an interaction term is present in a spa- 
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